Le gaz carbonique est un excellent gaz à effet de serre, comme 
on vous le serine tous les soirs au journal du soir-espoir, et au 
journal du matin-chagrin. 


Ce qui est dommage, c’est que personne ne vous dit selon 
quel mécanisme physique d’échange entre infra-rouges et 
molécule cet effet a lieu... 


Hi ! hi ! 


Comme le « gaz-carbonique-est-un gaz-à-effet-de-serre », il en 
résulte pourtant une amusante conséquence dont personne 
ne parle. 


Le Soleil nous envoie une quantité constante de rayonnement 
lumineux selon le spectre du corps noir. 


Donc, il nous envoie des infra-rouges, qui vous cuisent sur les 
plages. 


Mais comme le gaz carbonique est un gaz à effet de serre, -bis 
repetita placent-, il REPOUSSE donc la plus grande part du 
rayonnement infra-rouge d'origine solaire. 


Ainsi, le CO2 protége la Terre du rayonnement solaire. 


Hi! hi! hi! 


Gérôme XY Taillandier, Schweisthal. 


PS : et n’oubliez pas de vous cramer aux UV pour attraper votre 
mélanome cet été ! 


Dans ce qui suit, l'ai souvent utilise les termes 
apogee et perigee, au lieu de: aphelie et perihelie 
Vous corrigerez la ou C est necessaire. 


1 high eccentricity ellipsis 
In a case of HIGH eccentricity Earth trajectory, 


2 low eccentricity ellipsis Sun radiations are LOW on Earth 


EFFECTS OF EARTH ORBIT ECCENTRICITY ON SUN HEAT TRANSFER TO EARTH 


Since 200 OOO years at least, eccentricity has deemed, and is at present time very low. 
Heat transfer is then near its maximum. 


Eccentricity is by far the Milankovic parameter of outstanding importance 


Climate change is NOT due to carbondioxyde greenhouse effect, but to orbital heat transfer. GT 


Dans ce schéma la loi des aires de Kepler n'est pas respectée. 

La durée de séjour de la Terre dans le voisinage du Soleil n'est pas 
proportionnelle à la longueur du segment d'ellipse, 

mais diminue dans la région "chaude". 


vitesses de ks 
élevée 


1- La loi des aires montre que la vitesse de la Terre varie sur son orbite, 
plus élevée au périgée, plus faible à l'apogée. 


2 La quantité d'énergie émise par le Soleil est constante avec le temps 
3- La quantité d'énergie recue par la Terre varie en 1/r^2 avec la distance au Soleil. 


CONCLUSION: Plus l'excentricité de l'orbite de la Terre est faible, plus l'énergie reçue du Soleil est élevée 


Par conséquent, plus les températures de la Terre sont élevées; 
et pas le contraire ! GT 


Exentricité de l'orbite terrestre 


un Ы Periode intermédiaire 


climatic precession 
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Les variations de tous les paramètres. 
Voyez mon schéma. 


Le prétendu “réchauffement climatique" est dû à des gens qui 
ne savent pas raisonner sur la structure de l'excentricité de 
l'orbite de la Terre, et ne voient pas que, en raison de la loi de 
Kepler, la vitesse de la Terre prés du foyer solaire diminue le temps 
d'exposition au rayonnement, et que, de plus, la partie 
"éloignée" proche du second foyer a une exposition qui diminue 
en 1/r2, selon la loi bien connue de tous. 

Une excentricité FORTE DIIMINUE l'exposition de la Terre 

à la chaleur solaire. 

Géróme Taillandier, Schweisthal, 

2023 2 16 


La température sous abri, à Violay, lors de l'éclipse solaire partielle du 11/08/1999. 
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1 high eccentricity ellipsis 
In a case of HIGH eccentricity Earth trajectory, 


2 low eccentricity ellipsis Sun radiations are LOW on Earth 


EFFECTS OF EARTH ORBIT ECCENTRICITY ON SUN HEAT TRANSFER TO EARTH 


Since 200 OOO years at least, eccentricity has deemed, and is at present time very low. 
Heat transfer is then near its maximum. 


Eccentricity is by far the Milankovic parameter of outstanding importance 


Climate change is NOT due to carbondioxyde greenhouse effect, but to orbital heat transfer. GT 


RECHAUFF' MON PETIT CŒUR GELÉ ! 


Géróme Taillandier 


Selon une nouvelle trés répandue dans les journaux du matin 
et du soir, le gaz carbonique serait responsable d'un effet de 
serre qui provoquerait l'augmentation des températures 
atmosphériques. D'une part, cette affirmation confond 
chaleur et température. 


De plus elle n'explique rien. Si l'on admet, ce que je n'ai рй 
vérifier nulle part, qu'une molécule de gaz carbonique, CO2 
admet une charge de rayonnement infra-rouge, il est clair que, 
comme toutes les molécules, elle se désexcite au bout d'un 
temps de relaxation en émettant un rayonnement 
électromagnétique dans d'autres longueurs d'onde. Si l'on 
admet alors que, à un instant donné, il existe une proportion 
de CO2 dans l’atmosphère, si l'on exclut un apport nouveau, 
l'atmosphére se comporte comme un gaz en équilibre, dans 
lequel la chaleur des échanges internes est constante. 


Seul un apport de chaleur externe au systéme peut modifier 
cette température du systéme. 


*Considérons alors un modèle simplifié d'atmosphére, 
constituée d'une couche de gaz de 50 km de hauteur, et d'une 
surface de 510 millions de km2. 


*Négligeons les échanges avec le sol, en particulier avec le 
noyau terrestre, dont la chaleur doit pourtant bien se trouver 
quelque part en surface. 


*Les échanges à l'instant t dans la couche sont constants. 


*Un échange entre la couche et le vide entourant la Terre est 
envisagé. 


Pour un km2 de surface d'échange, la colonne de gaz de la 
couche est de 50 km3. 


*11 parait que la perte de chaleur et donc de rayonnement 
infra-rouge par échange avec le vide soit de 1 sur 50. 


SI l'atmosphère terrestre n'avait pas d'autres ressources, il est 
clair que l'espéce Sapiens serait congelée en trés peu de 
temps... 


Et ce n'est pas le CO2 qui la réchaufferait ! 


Par ailleurs, sous une éclipse totale de Soleil, on observe une 
perte de température de 3 *C environ durant prés d'une heure 
30, avec bien sûr une variation de cette perte ! 


Enfin, je rappelle que l'éruption volcanique islandaise de 1783 
a plongé les pays d'Europe du nord et de l'ouest dans un 
désastre alimentaire sans égal, qui durait encore en 1789. 


Ce n'est pas le CO2 du volcan qui les a réchauffés ! 


Je vous laisse le plaisir de consulter la quantité de chaleur 
recue du rayonnement solaire infra-rouge, qui vous permet de 
crever de soif dans un désert et de vous mettre à poil sur les 
plages pour pouvoir attraper votre mélanome, dû au 
rayonnement ultra-violet du Soleil. 


De deux choses l’une pour conclure. 


Ou bien le facteur de 1/50 proposé est suffisant pour 
accumuler la chaleur de l’activité terrestre, Sapiens ou non; 
et l’espèce Sapiens devrait être cuite depuis longtemps. 


Ou bien ce facteur est insuffisant; et les échanges 
électromagnétiques avec la molécule de CO2 ne sont en rien 
suffisants pour réchauffer le climat. 


En clair, le prétendu effet de serre des gaz est une foutaise, 
comme l'on déjà remarqué d'autres auteurs, qui soulignent 
que le CO2 atmo sur Terre était bien plus élevé durant les 
périodes glaciaires anciennes. 


Je vous donne rendez-vous pour patiner sur la Seine ou la 
Tamise dans quelques années, comme c'était le cas dans les 
années 1800. 


Géróme XY Taillandier, Schweisthal 


Saint Valentin 2023 
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Figure 4 - Tracé de la température de l'air près du sol 

(sous abri) mesurée toutes les minutes entre 9 h et 15 h UTC 
le 11 août 1999, pour trois stations de Météo-France 

situées dans la bande de totalité. 

Le trait vertical discontinu gras indique l'heure centrale 

de l'éclipse totale. 

Les traits verticaux discontinus minces 

indiquent l'heure du premier contact (début de l'édipse partielle) 
et du demier contact (fin de l'éclipse partielle). 

La ligne en pointillés et la flèche montrent comment on a estimé 
le déficit de température dû à l'éclipse. 

a - Station météorologique de Saint-Quentin (Aisne) ; 


b - Station météorologique de Beauvais (Oise) ; 


€ - Station météorologique d'Abbeville (Somme). 


La température sous abri, à Violay, lors de l'éclipse solaire partielle du 11/08/1999. 
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THE SUN AS А LOW-FREQUENCY HARMONIC OSCILLATOR 
P. E. DAMON and J. L. JIRIKOWIC 


Laboratory of Isotope Geochemistry and The NSF-Arizona Accelerator Facility for Radioisotope 
Analysis, Department of Geosciences, The University of Arizona, Tucson, Arizona 85721 USA 


ABSTRACT. Solar activity, as expressed by interplanetary solar wind magnetic field fluctuations, modulates the 
atmospheric production of “C. Variations of atmospheric “C can be precisely established from the cellulose within annual 
tree rings, an independently dated conservative archive of atmospheric carbon isotopes. AC time series interpretation shows 
that solar activity has varied with a recurrence period of 2115 * 15 (9596 confidence) yr (Hallstattzeit) (Damon & Sonett 
1991) over the past 7160 yr. From a non-stationary oscillation solar activity hypothesis, 52 possible spectral harmonics may 
result from this period. Damon and Sonett (1991) identify powerful harmonics such as the 211.5-yr (Suess) and the 88.1-yr 
(Gleissberg) cycles as independent fundamental periods. These stronger harmonics appear to modulate the 11-yr (Schwabe) 
sunspot cycle. Variations in the solar magnetic field, thus, may respond to longer period variations of the solar diameter 
envelope (Ribes et al. 1989). Such variations would affect solar radiative energy output and, consequently, change total solar 
irradiance (Sofia 1984). 


HARMONIC ANALYSIS OF THE “C SPECTRUM 


Spectral analysis of the available high-precision AC data by Stuiver and Braziunas (1989), using 
the maximum entropy method (MEM), led them to postulate a fundamental 420-yr oscillation 
related to changes in the solar convective zone. After increasing the model autoregressive order 
(AR, equivalent to the length of the prediction error filter, for the purposes of this paper), Stuiver 
and Braziunas observed that additional periods became important. They note that the 420-yr period 
splits into periods of 504, 355 and 299 yr. We have observed that this apparent splitting of the 420- 
yr peak results from the limited frequency resolution, and to oversmoothing by the MEM spectral 
method when too small an AR order is chosen. Our experience suggests that frequency resolution 
for MEM (AR = 160) spectral analysis lies between 1/x and 1/1.25т, where т is the length of the 
time series. Although the analysis above seems to suggest subjectivity, MEM analyses provide 
excellent frequency resolution when comparing between many AR orders and with the results of 
alternative spectral estimate methods (such as the discrete Fourier transform (DFT), multi-taper 
estimates, maximum-likelihood-Bayesian methods and singularity spectral-analysis estimates). By 
varying the number of sampled points and AR order, it becomes apparent that the 420-yr peak does 
not split, but each of these peaks is always present with relative peak heights varying with chosen 
AR order. 


However, these low-frequency variations remain poorly resolved, whatever the spectral method 
employed, until either the data time series lengthens, or a hypothesis develops relating the poorly 
resolved low-frequency variations to better-resolved higher-frequency variations (see Damon & 
Jirikowic 1992, for details). We have used the simplest non-stationary model, a harmonic oscillator, 
as our working hypothesis. The well-resolved higher-frequency variations represent a set of 
harmonic and combination overtones of lower-frequency variations. For example, concerning the 
~504-yr period, its 2nd-13th apparent overtones observed in spectral analyses of the last 7160 yr 
of the high-precision bidecadal AC calibration time series (Stuiver & Kra 1986) establish the 
period more precisely as 526 yr. Similarly, the 9 apparent overtones of the 420-yr peak define it 
as 424 yr; the 8 apparent overtones of the 355-yr period establish it as 356 yr, and analysis of the 
299-yr period with 6 apparent overtones established it more precisely at 300 yr. Further, each of 
these longer periods appears to modulate shorter periods producing peak triplets, quintets and 
higher. Further analyses show these spectral peak frequencies may be related as elements of the 
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harmonic series T/n, T/2n, T/3n, . . ., where n is any positive integer, 1, 2, 3,... Assuming n = 
1, and finding an approximate least-common multiple to obtain the shortest possible estimated 
period of the fundamental, T, we obtained (keeping 4 significant figures to reduce truncation 


errors): 
526 x 4 = 2104 
424 x 5 = 2120 
356 x 6 = 2136 
300 x 7 = 2100 


T = 2115 + 9 yr (o) 


Thus, the supposed fundamentals may be harmonic overtones of a ~2115-yr Hallstattzeit period 
(Damon & Sonett 1991). By hypothesis that the Hallstattzeit period may have 52 harmonic 
overtones present in the spectral analyses (limited only by the 40-yr Nyquist period of the 
bidecadal AC time series), the fundamental period may be established from the observed peak 
power frequencies with a precision of 0.1 + 3.1 yr (95% confidence level). The agreement of the 
Hallstattzeit overtones with reference to both the DFT and MEM (AR = 160) spectra (Fig. 1) 
strongly suggests that the causal phenomenon is independent of the spectral method. Thus, 
modulation side-lobe “leakage” due to spectral window truncation does not seem to influence 
greatly our results. 


PROBABILISTIC VALIDATION OF THE HALLSTATTZEIT OVERTONES 


The probabilistic! significance of this agreement may be investigated in several ways. For 
example, of the 20 DFT frequencies with spectral power > 95% x?-confidence, 7 have frequencies 
matching the hypothesized Hallstattzeit harmonic overtone frequencies within the frequency 
resolution. From the discrete, multinomial distribution of these sets of frequencies, the probability 
of equaling the performance of the Hallstattzeit harmonic hypothesis with 52 randomly chosen 
frequencies is 0.12%. The likelihood of the harmonic hypothesis reaches a maximum near a 
fundamental period of 2115 yr, and remains above 95% from 2120 to 2105 yr. From Bayesian 
analysis, we are able to infer from the null hypothesis a priori probability distribution of random 
matching that the a posteriori probability of random matching remains zero within 95% confidence 
for hypothetical fundamental periods between 2130 and 2100 yr. The final test employed the 
second order norm between the observed spectral peak frequencies and the hypothesized overtone 
frequencies. A Hallstattzeit fundamental period of 2116-yr minimized the second order norm with 
a 2o (~95%) range from 2128 to 2104 yr. Thus, the solar harmonic oscillator hypothesis with a 
fundamental Hallstattzeit period between 2130 and 2100 yr provides a compelling explanation of 
much of the observed AC time series spectrum. 


MODULATION OF THE SOLAR SCHWABE CYCLE BY THE SUESS AND GLEISSBERG HARMONICS 


A harmonic oscillator would be a particular solution to the general differential equation 


DO + f[D(t)] + g[D()] = F() (1) 


Because the hypothesized 2115-yr Hallstattzeit variation recurs only three times in the МС calibration time series, and the 
harmonic overtones may not have independent statistical distributions, non-parametric probability analyses may be more 
reliable and informative than traditional parametric statistics. 

2At any given frequency, the MEM absolute power depends upon the chosen AR order. Thus, conventional confidence-level 
calculation analysis of MEM spectral peaks is inappropriate without independent a priori determination of the AR order. 
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Fig. 1. A. DFT power spectrum of the high-precision А!“С data. Vertical lines represent the Hallstattzeit fundamental and 
52 harmonic overtones. The heavy curve estimates the lower 95% x’-confidence limit and the “red-noise continuum" likely 
due to non-stationarity and time persistence. B. MEM (AR = 160) spectrum. Note that although much of the spectrum 
trends as a red-noise continuum, 20 spectral peaks extend above the 95% X-confidence limit, suggesting substantial 
periodicity. C. compares the observed and predicted spectral peaks. The error bars on the top of the observed lines represent 
the frequency resolution. 


where D(t) represents solar activity, f, a damping function, g, a restoring function and F(t), an 
external forcing. The harmonic oscillator particular solution of Equation (1) may take the form of 
a Fourier sum of harmonics 


Dit) = A, + У)" [Acos(ot-b)-Bsin(wt-)] (2) 


where œ; are harmonic frequencies, A, and B, are amplitude coefficients and $, are phase 
coefficients. Note that A, B; and $; may be complex functions of both the harmonic frequencies 
and time. 


Certain harmonics of the Hallstattzeit with observed power greater than expected for a simple 
system of harmonics act as apparent fundamentals (i.e., 211.5-yr Suess and 88.1-yr Gleissberg 
variations). These periods appear in the DFT and MEM (AR = 160) spectra of the annual Wolf 
sunspot indices as harmonic overtones and modulating periods (Fig. 2). For example, the first two 
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Fig. 2. А. РЕТ power spectrum of the annual Wolf sunspot index. Vertical lines represent the 211.5, 88.3 and 11.1-yr 
periods, their overtones and sidebands produced by amplitude modulation of the 11.1-yr carrier period (after Sonett 1982). 
The heavy curve estimates the lower 95% x?-confidence level and the red-noise continuum. B. MEM (AR = 160) spectrum. 
Vertical lines represent observed peaks. Modulation sidebands appear in the annual Wolf sunspot index spectra, features 
not seen in the AC spectra in Figure 1. C. compares the observed spectral peaks with those predicted from the amplitude 
modulation model of the 11-yr cycle. The error bars on the top of observed lines represent the frequency resolution. 


harmonic overtones of the 211.5-yr Suess period cannot be resolved from the 88.1-yr Gleissberg 
period, and not all side bands of the 11.1-yr Schwabe cycle carrier modulated by these longer 
periods and their overtones can be resolved from one another (Fig. 2). Higher-frequency overtones 
appear, including four overtones of the 11.1-yr carrier (5.55, 3.70, 2.78 and 2.22-yr, see Fig. 2). 
The 211.5-yr period and its 2nd (105.8-yr) and 4th (52.9-yr) overtones, combined with the Gleiss- 
berg period (88.1-yr), appear to modulate strongly the 11.1-yr carrier. A reasonable fit to the 
observed annual Wolf sunspot indices can be obtained, using a squared amplitude modulation 
model (after Sonett 1982) by allowing these periods (211.5, 105.8, 52.9 and 88.1-yr) to amplitude 
modulate the 11.1-yr carrier signal (see Fig. 3, Table 1) 


к-У; 2 a,cos(wt + ф,)о соз (0 +ф,) . (3) 
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Fig. 3. (-—-) simulation of the (—) annual Wolf sunspot indices by modulation of the 11.1-yr carrier by the 211.5, 105.75, 
88.125 and 52.875-yr periods with the amplitudes and phases shown in Table 1. Note that the 211.5-yr periodicity has a 
low modulation coefficient. Also note that the model underestimates sunspot indices during the late 19th century solar 
activity minimum. Because the Hallstattzeit gate strongly affects the 211.5-yr harmonic and this millennium's Hallstattzeit- 
gating epoch ends before the late 19th century, the weak 211.5-yr and the weak late 19th century solar minimum would 
support the Hallstattzeit non-stationarity hypothesis (see text for details). 


TABLE 1. Model Parameters for Modulation of the 11.1-Year Schwabe Carrier 


Index Period a $ 
1 52.875 0.66 0.62 
2 88.125 0.43 -1.82 
3 105.75 1.31 1.10 
4 212.5 0.01 1.33 
c 11.1 65.4 1.46 


Note that the cosine series modulation function expressed in Equation 3 is a member of the family 
of harmonic oscillator solutions given in Equation 2. Notably, the model fits the late 19th century 
minimum poorly. Mean sunspot indices do not fall to the very low levels predicted by our model. 
Also note that the modulation coefficient, alpha, for the 211.5-yr period is quite small. Both events 
may be explained by the non-stationarity present in the Wolf sunspot index data set. The 
Hallstattzeit epoch encompassing the Wolf, Spórer, Maunder and Dalton solar activity minima 
wanes in the first half of the sunspot index data. Thus, the stationary model assumption has been 
compromised. The Hallstattzeit variation also appears to modulate strongly the 211.5-yr period 
(Damon & Sonett, 1991). A more sophisticated model would include the Hallstattzeit gating. 


HALLSTATTZEIT GATING: PHYSICAL MANIFESTATION AND PALEOCLIMATIC CONSEQUENCES? 


The Hallstattzeit variation appears.to act as a gate (Damon 1988) allowing more intense century- 
scale oscillation during the high AC phase of the cycle (e.g., the Maunder and Hallstattzeit low 
solar activity epochs; Schmidt & Gruhle 1988). If we assume the wave form for the Hallstattzeit 
period, shown in Figure 4, with the high AC phases appearing as the square-wave gate portions, 
the square wave generates the harmonic overtones and strongly modulates the 211.5, 88.1-yr 
periods. Figure 4 presents the AC variations after the removal of a very long-term trend ascribed 
to changes in the Earth's magnetic dipole moment (see Damon & Sonett 1991 for review). The 
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Fig. 4. (---) wave-form model generating the Hallstattzeit harmonics and amplitude gating of the (—) century-scale AC 
variations compared to the detrended high-precision AC time series. Note that the weaker century-scale oscillations appear 
modulated by variations of the Earth’s magnetic field as described in the text, yet the stronger gated variations appear 
comparatively unmodulated. 


trend can be closely fit by a sine curve with a 10,780-yr period for the last eight millennia, but 
deviates markedly during the early Holocene (Damon, Cheng & Linick 1989). The Earth’s 
magnetic dipole moment is weak during the negative portion of this cycle before 2600 BC and 
stronger during the positive part of the cycle between 2600 BC and AD 1650. The effect of this 
geomagnetic field variation can be seen by relative suppression or enhancement of the century-scale 
AC variations. This amplitude modulation of AC variations by the geomagnetic field appears 
weak during the high AC gating portion of the Hallstattzeit cycle. The Hallstattzeit variation may 
have also affected paleoclimate by changing solar output and the Suess and Gleissberg cycles may 
have contributed in part to the global warming trend evidence since the last decade of the 19th 
century (Damon 1989; Damon & Jirikowic 1992; see Damon & Sonett 1991 for review). 


The physical manifestation of the Hallstattzeit solar variation has yet to be conclusively determined. 
Ribes (1990: 96), reviewing the evidence for periodicities in the solar diameter, concludes that, “the 
Sun’s apparent radius was larger in the deep of the Maunder Minimum (late in the Little Ice Age, 
when practically no sunspots were seen) than it was at the end when solar activity has resumed 
(from 1705 onwards). An expansion of solar envelope of about 3 arc seconds on the diameter 
would be consistent with cooling of the envelope and with a slowing down of surface rotation.” 
With relation to Ribes’ (1990) interpretation, the Hallstattzeit wave form in Figure 4 suggests that 
the Sun slowly expands, lingers at maximum diameter and then contracts. This structural modu- 
lation partitions potential energy and radiant energy (Sofia 1984) between two intransitive states. 


CONCLUDING REMARKS 


The harmonic overtone periods may represent a dynamic-system response and other non-stationary 
phenomena or a mathematical manifestation of the non-sinusoidal Hallstattzeit variation. The 
observation of variations, such as the =88-yr Gleissberg cycle, strongly suggests that several 
variations also have a physical reality beyond being harmonic overtones of the fundamental 
Hallstattzeit variation. Such secular solar variations would have obvious implications as possible 
factors in global climate change. Associating these Hallstattzeit solar variations with climate 
variations may require a time-variant climate model. Further, the apparent non-stationarity of the 
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solar variations prohibits extrapolating linear regression models calibrated over a single 11-yr cycle 
between solar activity and solar radiant output beyond a Hallstattzeit variation. Predictions based 
solely upon the cyclicity of A'*C variations, without consideration of these Hallstattzeit non- 
stationarities, may be more misleading than informative. Understanding both the modeling and 
geophysical implications of the Hallstattzeit variation would advance our ability to predict solar 
activity and understand the Sun as a harmonic oscillator. 
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Solar Forcing of Global Climate Change 
By PAUL E. DAMON AND JOHN L. JIRIKOWIC 


Department of Geosciences, University of Arizona, Tucson, AZ 87521, USA 


Using present global warming and paleoclimatic records from climatically sensitive regions as 
a frame of reference, we infer that global temperature changes did not exceed +0.5°C during 
the current interglacial or +2.0°C during the last glacial period. In order to completely explain 
such fluctuations by solar irradiance changes, all solar and terrestrial factors must be optimized. 
Since this is unlikely, we conclude that solar forcing of pre-anthropogenic climate change is a 
significant and perhaps dominant factor but other processes must also be significant. Solar 
irradiance changes alone cannot result in the temperature decrease required to change from 
interglacial to glacial as suggested by Opik (1965). 


1. Introduction 
1.1. Glacial and interglacial climate 


During the Pleistocene (1.8 Ma to present), the Earth's climate has varied between two 
states, glacial and interglacial. This can be seen in the oxygén isotope stratigraphy 
in Figure 1. The meteorological cycle has been likened to a multiple-stage still with 
6 to 7 theoretical distillation plates (Epstein 1958). The amount of fractionation in 
the hemispheric distillation process is a function of altitude, latitude and temperature 
gradient from equatorial to polar regions. The variations of the 180/160 ratio during 
the hydrologic cycle are measured in %o relative to mean seawater. With reference to 
the Byrd core, Antarctica, it can be seen that 6180 fluctuates around —40% from about 
110 ka to 18 ka, gradually increasing after 18 ka to fluctuate around a new steady-state, 
—34%o, from about 10 ka to present. The increase in 8180 from —40%o to ~34%o is the 
result of global warming and coincident with rising sea-level as the continental ice masses 
melted (Fairbanks 1989). The Camp Century, Greenland core shows 6:80 variation 
around the same two steady states but during the transition form 18 ka to 10 ka, the 
warming is temporarily interrupted by a return to glacial conditions during the Younger 
Dryas from about 12,900 BP to 11,600 BP. Thus, the hydrologic cycle is responding 
to both regional and global climate changes. The North Atlantic is more sensitive to 
climate change during the glacial period as can be seen by the larger fluctuations of 
6180 in Greenland ice cores when compared to Antarctic ice cores but the two states are 
global in nature. In general, global climate during the Pleistocene can be modeled as a 
set of systems fluctuating around the glacial and interglacial steady states separated by 
transitional episodes. 6180 fluctuations around steady state are three to four times less 
during the interglacial compared to the glacial epoch. 


1.2. Regional vs. global climate change 


According to estimates from the CLIMAP project (CLIMAP Project Members, 1981), 
the Laurentides ice mass reached a thickness of 3800 m and the northwestern European 
ice mass reached a thickness of 2500 m, whereas the Greenland ice cap reached an 
intermediate thickness of about 3200 m. Using an empirical equation from Dansgaard 
(1964), we can estimate the glacial age temperature at Camp Century to have been 
& —38°C. The Northern Hemisphere continental areas covered by ice must have been 
at similar temperatures, warming to their present annual temperatures. As a robust 
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FIGURE 1. Ice core 5180 (Camp Century, Greenland, after Dansgaard et al. 1982 (1a); Byrd 
Station, Antarctica, after Johnsen et al. 1972, (1b)). Note that both isotopic records show . 
distinct modes between glacial and interglacial states. For the glacial state (20 ka-13 ka), 
6180 varies by an amplitude of 2%o on a century time-scale about an average value of zz —40%o. 
During the interglacial (9 ka-present), 61° О varies by an amplitude of «0.5960 around an average 
of © —29%o for Camp Century and zz —33%o for Byrd Station. The decreased 5:80 variability 
suggests a decrease in sensitivity to climate forcings. Camp Century shows greater glacial to 
interglacial contrast due the influence of surface ice and snow expanse during the glacial period. 
A sharp increase in the Camp Century record contrasts with the more gradual increase in 
the Byrd Station record. This may reflect the ameliorating influence of the Southern Ocean. 
The pronounced Younger Dryas negative deflection (2210.5 ka) present in the Camp Century 
may correlate with a much smaller feature (7:10 ka) in the Byrd Station record (calibrated 
ages are about 1700 years older as in the text). The Younger Dryas event may have affected 
North Atlantic circulation greatest and thus adjacent continental indicators have more profound 
signatures. Thus, 5180 records contain both global and region climate signals. 
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Region AT (°С) Area of Globe (%) Contribution to Global 
Mean weighted by area (°C) 

Presently Glaciated —10 3 —0.3 

Not glaciated at present —40 5.7 —2.3 

Other continental areas —6 20.5 —1.2 

The world oceans —2.3 70.8 —1.6 

Hansen et al., 1984 У = —5.4 

СОМ Model: У = —5.3 


TABLE 1. Temperature changes (AT) from present to 18 ka. 


but not rigorous estimate of the temperature increase for those areas, we suggest about 
40°C. For the areas still covered by polar ice, we estimate a change of about 10°C. From 
a paleoclimatology summary by Crowley & North (1991), we take the average change at 
mid-latitudes to be between 5-8? C and 4-5°C for tropical continent and areas under the 

' influence of maritime air masses. Let us assume an average of ~6°C for continental areas 
not covered by ice sheets. CLIMAP Project (McIntyre et al. 1976) estimated an average 
change for the global oceans of only 2.3°C but with great variation between different 
ocean basins. For example, in the Northern Hemisphere, they estimated temperature 
changes varied from 0.8°C for the Indian Ocean to 3.8°C for the North Atlantic. We 
summarize and weight these changes by areas in Table 1, arriving at a global change 
of —5.4?C from the glacial to the present state. This compares favorably with more 
rigorous GCM (global circulation model) modeling by Hansen et al. (1984) and serves 
to ilustrate several pertinent points. First, the global oceans serve as a great source 
and sink of heat moderating climate. Second, 8.796 of the area of the Earth accounts 
for &50% of the global warming during the transition from glacial to interglacial states. 
Third, there were large regional differences in temperature changes around the global 
mean. 

As previously noted, oxygen isotope stratigraphy demonstrates that the hemispheric 
temperature fluctuations during the past 10,000 years (the Holocene) are relatively small. 
The linear empirical relationship between 6195O and temperature (Bradley 1985) suggests 
that these fluctuations do not exceed +0.6°C (or 1.2? C peak-to-peak). If this is correct, 
global warming during the past 100 years of 0.5? С based on instrumental records (Folland 
et al. 1990) should be a very significant event relative to previous events during the 
Holocene. Of course, regional temperature variations of sub-hemispheric scale have been 
both lesser and greater. For example, there was an insignificant increase in temperature 
during the last century for the contiguous United States comprising roughly 1.696 of 
global area (Hanson et al. 1989). However, from data in Trends, 1991 (Boden et al. 
1991), we calculate an increase of (1.2+0.2)°C(s) from solar cycle 14 (1901-1913) to 
solar cycle 21 (1976-1987) for the southern deserts which is one of 23 regions of the 
U.S.A. reported in Trends, 1991. These sorts of regional variations of climate change 
together with conceptual and dating problems make it difficult to evaluate past changes 
in Holocene global temperatures. The instrumental record coupled with GCM studies 
can be helpful in resolving this problem. Specifically, in order to evaluate the feasibility 
of solar forcing of global climate, we need to know the magnitude of global Holocene 
temperature changes to be forced. As noted above, oxygen isotope studies of polar ice 
suggest global changes equal to or less than +0.6°C during the Holocene. This suggestion 
needs further justification. 

GCM studies suggest that global temperature change will be greatest at high altitudes 
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and high continental latitudes (Wang et al. 1991). High altitude temperatures are 
governed by evaporation of surface ocean water with transfer of 540 calories per gram 
of latent heat to higher latitudes. Albedo changes, as the result of changes in the ice 
and snow cover, produce high-latitude continental sensitivity to global climate forcing. 
Evidence for climatic events like the Medieval Warm Epoch and Little Ice Age come 
from North Atlantic border-lands and Europe above 45?N. It has been possible recently to 
demonstrate by study of montane glaciers that these climatic events are global rather than 
merely regional in scope (Rothlisberger 1986). Montane glaciers also provide perspective 
concerning the relative intensity of 20th Century global warming. 

Hasenrath & Kruss (1992) have compared glacial maps for Mount Kenya (elevation 
4.6 to 5 km). They estimate that the Mount Kenya Glacier retreated by 40% in area 
from 1963 to 1987. The retreat of the terminus of the Speke Glacier in the Ruwenzori 
Range of Uganda has accelerated from an average of 2:2 m/yr between 1958 and 1977 to 
more than 11 m/yr between 1977 and 1990 (Kaser & Noggler 1991). The areal extent of 
worldwide montane glaciation has been reduced by about 5096 during the last century 
(Thompson 1991). Thompson et al. (1989) have studied the Dunde Ice Cap on the 
Qinghai-Tibetan Plateau at a summit elevation of 5235 m. According to these authors, 
in that region "the last 60 years were apparently one of the warmest periods in the entire 
record, equaling levels of the Holocene maximum between 6000 and 8000 years ago" (p. 
474). 

Evidence for high-latitude warming comes from dendroclimatologic studies and mea- 
surements of the temperature gradients in bore holes. Innes (1991) has documented high- 
altitude and high-latitude increases in 20th Century tree growth from about 30 studies of 
trees from Northern Europe and North America. In contrast to Thompson's conclusion, 
Innes states, “It appears that neither the northern extent nor the present growth rates of 
boreal forests have reached the extent of rates characteristic of the mid-Holocene, when 
the climate was demonstrably warmer." However, from their dendroclimatologic study 
of high-latitude boreal sites in Canada and Alaska, Jacoby & Arrigo (1989) conclude 
that "the global warming trend over the past century of increasing СО» has exceeded 
the recent (past three centuries) level of natural variability (pre-anthropogenic) of the 
climate system". In a study of old-growth Huon pines from a sub-alpine site in western 
Tasmania, Cook et al. (1992) concluded that, *Although the most recent 25-year period 
is warmer than any comparable period over the past 1090 years, it is not yet statisti- 
cally unprecedented" (p. 215). Graybill & Shiyatov (1992) have reconstructed June-July 
temperatures for ad 961-1969 using tree-ring width variation of Siberian Larch from the . 
Polar Urals (65-67?N). They concluded that “the relatively sustained rise in tempera- 
ture from the 1880s to the 1960s has the greatest magnitude in the record for such an 
extended period (1.6°C difference between the 20-year intervals centered on 1870 and 
1960)" (p. 411). Note that for this climate sensitive area, this recent regional warm- 
ing event is threefold greater than 20th Century global warming of ~0.5°C. Bore hole 
temperature gradients from Central and Eastern Canada (Beltrami et al. 1992; Shen & 
Beck 1992; Wang et al. 1992) as well as from the Alaskan Arctic (Lachenbruch & Vaughn 
1986; Kakuta 1992) confirm the temperature rise from the late 19th Century through 
the late 20th Century. The temperature rise in Canada was typically 1-3°C and greater 
at Prudhoe Bay, Alaska where Kakuta (1992) estimates warming by 4°C from ad 1853 
to ad 1973. 

Summarizing, the modern instrumental record has made it possible to estimate 20th 
Century global warming of about 0.5°C (Folland et al. 1990). There is insufficient data 
to establish directly global climate change prior to the 19th Century. The existing data is 
biased toward relatively high-latitude, climatically-sensitive regions of the North Atlantic. 
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However, our comparison of the recent warming event with the paleoclimate record within 
climatically-sensitive regions suggests that this recent event exceeds in intensity natural 
global climate variability during the last millennium and is of comparable magnitude to 
the Holocene maximum between 6000 and 8000 years ago. Oxygen isotope stratigraphy 
also demonstrates the relative global stability of pre-anthropogenic climate. Solar forcing 
needs only account for global temperature variations of no more than +0.5°C to account 
for all the natural climate variability during the last millennium. The Holocene maximum 
may exceed a global temperature increase of 40.5? C. However, our previous estimate of 
+0.6°C variability from the polar oxygen isotope record is probably a conservative upper 
limit for global temperature change during the Holocene. 


2. Solar forcing of Holocene climate 
2.1. The Sun as a low-frequency harmonic oscillator 


The production of cosmogenic 14C is modulated by solar activity via the magnetic field 

associated with the solar wind. The relationship is inverse and, so, the various histori- 
cally observed solar minima such as the Maunder Minimum occur as !*C maxima and 
the Medieval Solar Maximum from ad 1100-1250 occurs as а 1*C minimum in the spec- 
trum of !*C variations. The reader may refer to two papers for recent reviews of the 
subject (Damon & Sonett 1991; Stuiver et al. 1991). Damon & Sonett (1991) concluded 
that medium- and short-term variations with periods of 2:208 (Suess), 88 (Gleissberg), 
and 11 (Schwabe) years are in large part produced by solar wind modulation. A longer 
period of &2300 years (Hallstattzeit) seemed enigmatic but appears to 100% amplitude 
modulate the Suess cycle. The &2300-year period seems to be associated with recurrent 
Little Ice Ages (Schmidt & Gruhle 1988). Stuiver & Braziunas (1991), using the maxi- 
mum entropy spectral method (MEM), chose a 420-year period as a fundamental with 
harmonics related to changes in the solar convective zone and noted that it splits into 
three periods at higher AR order. We note that the apparent splitting is an artifact of 
the MEM method. The peaks are always present but vary greatly with the autoregres- 
sive (AR) order. Too low an AR order results in oversmoothing. We deduced from both 
DFT and MEM power spectral analyses that the four periods (526, 424, 356 and 300 
years) are the 4, 5, 6, 7 harmonics of the poorly resolved Hallstattzeit period that we 
estimated from its harmonics to be 2120 years (Damon & Jirikowic 1992a,b). Further 
analysis yielded 212 years for the Suess period with strong (29596 confidence) harmonics 
at 106 years (2nd harmonic) and 53 years (4th harmonic). 

The harmonics of the Suess period combined with the Gleissberg period appear to 
strongly modulate the 11.1-year (Schwabe) cycle. Furthermore, a reasonable fit to the 
annual Wolf Sunspot Indices can be obtained using a squared amplitude modulation 
model (after Sonett 1982) by allowing these periods (212, 106, 53, 88.3 years) to amplitude 
modulate the 11.1-year carrier (Damon & Jirikowic 1992a). 

Subsequent to the work summarized above, one of the authors (Jirikowic) added Thom- 
son's (1990) multi-taper power spectral analysis technique to the techniques made avail- 
able to us by Prof. Charles Sonett here at the University of Arizona (periodogram, 
Bretthorst-Bayesian, Maximum Entropy [MEM], Discrete Fourier Transform [DFT]). 
The analysis of the latest version of the high-precision 14C data set (Stuiver & Reimer 
1993) is shown in Figure 2. The Hallstattzeit fundamental and first nine overtones (2- 
10fx) clearly resolve at 29576 confidence (harmonic F-test). The Suess period (fs and 
10fg) has unusually high power spectral density and appears to act as a fundamental 
as well as à harmonic of the Hallstattzeit. The accounts for the high spectral power 
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FIGURE 2. Power spectral density estimates using Thomson's Multi-taper Method (Thom- 
son 1990) and three autoregressive modeling orders: 120 (dotted), 160 (solid) and 200 (dashed) 
from the 1993 Radiocarbon Calibration tree-ring 5'*C data (Stuiver & Reimer 1993). Annotated 
peaks show assignment to hypothetical harmonics to the fundamental Hallstattzeit (subscript 
H), Suess (subscript S) and Gleissberg (subscript G) variations. All these peaks have greater 
than 9596 F-test confidence. Although the Hallstattzeit is poorly resolved, the annotated har- 
monics combine to suggest a fundamental period of 2108 + 15 years. Earlier spectral analyses 
(Damon & Jirikowic 1992b) showed a strong fourth harmonic of the Suess that is not apparent 
here. 


spectral density of its second harmonic (2fs, >95% confidence). However, its fourth har- 
monic appears as a relatively weak doublet, unlike the MEM and DFT techniques (>95% . 
confidence). The Gleissberg period is prominent at 29576 confidence (88.1 years). 

Neglecting the poorly resolved second Hallstattzeit harmonic but using harmonics 3-10 
plus 14, 15, 17, 18 (all with 29596 confidence) to calculate the Hallstattzeit fundamental, 
we obtain 2108+15 (2 years compared with 2115+18 (д) years previously obtained using 
the combined MEM and DFT methods (Damon & Jirikowic 1992b). 

The enigmatic Hallstattzeit period may be divided into two episodes. For 221520 
years, it is associated with high solar activity and weak periodicity followed by a ~600- 
year episode in which solar activity is low on the average and the Suess period and 
its harmonics as well as the Gleissberg period are strong. These 600-year episodes are 
associated with Little Ice Ages. Thus, when the Medieval Solar Maximum ended at 
ad 1250, the Little Ice Age began accompanied by a series of profound solar minima 
and terminated after about 600 years toward the middle of the 19th Century when the 
Modern Solar Maximum began (Jirikowic & Damon 1994). 

Instrumental analysis has confirmed a significant change in solar irradiance of 0.1% 
during the 21st solar cycle (Willson & Hudson 1988, 1991). Paleoclimate data suggests 
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Change in А Forcing AfTea (°C) ATea (°C) ATea (°С) 


Solar Solar! (tropopause  Glacial? Present Glacial? Future? 
Condition Irradiance level) 2.4°C 1.0°C 0.6°C 
Present Sun 0.0 0.0 0.0 0.0 0.0 
AD 1980-1986 

No plage, spots —1.5 —0.36 —0.9 —0.4 —0.2 
or network 

No-cycling stars —2.4 —0.58 —1.4 —0.6 —0.3 
(average) 

Minimum possible —3.5 —0.84 —2.0 —0.8 —0.5 
Sun 


"From Lean et al. (1992). 
2Calculated from Hansen et al.(1984). 
*Calculated from Hansen et al.(1984) assuming variations are negligible on an ice-free Earth. 


TABLE 2. Decreasing equilibrium temperature (T.,) as a function of solar condition. 


that longer periods such as the Hallstattzeit, Suess and Gleissberg periods may also be 
associated with changes in solar irradiance (e.g., see Eddy 1977). Further evidence comes 
from solar and astrophysics. 


2.2. The Sun as a variable star 


As discussed elsewhere in this colloquium's proceedings (see Lockwood 1994), 74 solar- 
type stars have been observed at the Mount Wilson Observatory. Of these 74, 13 have 
been studied monthly since 1966 and nightly since 1980 (Baliunas & Jastrow, 1990). Only 
sporadic measurements were made since 1978 for the other 61 stars. Nine of the 13 stars 
were observed to show a range of magnetic activity similar to the solar cycle whereas four 
that showed no variation of magnetic activity with time were assumed to be analogous 
to the Sun during the Maunder and other solar minimum. The nine magnetically active 
stars had regular activity cycles like the Sun. Their brightness varied by up to 0.496 
compared to 0.196 for the Sun during cycle 21. Baliunas & Jastrow (1990) state, “Our 
results suggest that when a solar-type star enters or leaves a Maunder Minimum, its 
brightness may change by more that the 0.1% observed during the Sun's 11-year cycle” 
(p. 522). Lockwood et al. (1992) presented a compilation of eight years of observations 
of 33 stars at the Lowell Observatory. They report that the stars whose mean activities 
were closest to the Sun's varied in brightness by 2.5 to 10 times the variation in brightness 
of the Sun. They state, "This raises the question of whether the Sun's cycle irradiance 
variation is always so small or whether larger changes occur in other cycles or over longer 
time intervals" (p. 655). 


2.3. Solar variability and the Maunder Minimum 


Based on current knowledge of the Sun and Sun-like stars, Lean et al. (1992) estimated 
the change in solar irradiance and forcing at the top of the atmosphere for various 
conditions of the Sun (see Table 2). Assuming, as did Baliunas & Jastrow (1990), that the 
Maunder Minimum was equivalent in change of irradiance to the average solar-like non- 
cycling stars yields a 2.4 change in irradiance. Also, assuming that a global equilibrium 
temperature decrease of 1?C was required to explain Little Ice Age cooling, they conclude, 
*.. that solar forcing may have contributed some, but not all, of the observed warming 
from the Little Ice Age to the present time" (p. 1594). 

Keeping in mind that the &0.5°C rise in temperature in recent time began from a 
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Feedback Gain AT (°C) 
Unity 0 -0.6 
Water vapor-clouds 0.415 —2.2 
Land ice 0.169 —0.9 
Sea ice 0.132 —0.7 
Vegetation 0.057 —0.3 
СО» 0.113 —0.6 


gi = 0.886 Y ATA = —5.3 


TABLE 3. Gains (gi) and temperature increments resulting from feedbacks at the Ice Age 
Maximum (18 ka). 


relatively cool global temperature in the late 19th Century, and this rise approaches or 
exceeds previous such warming events during the Holocene (see 1.2 above), we suggest 
that the estimate of 1°C decrease in global cooling exceeds what is required to explain 
the Little Ice Age cooling during the Maunder Minimum. We concluded that less than 
+0.5°C global cooling is required to explain climate variability during the last millen- 
nium. The global temperature rise from the Maunder Minimum to ad 1980-1986, based 
on the paleoclimatic record and instrumental measurements discussed, is greater than 
0.5? C but less than 1°C. 

In order to evaluate the temperature decrease resulting from a given change in solar 
irradiance (Table 2), various positive feedbacks must be considered. For our analyses, 
we refer to Hansen et al. (1984). Their paper included analysis of the various gains (g;) 
and temperature increments (AT;) including those for water vapor-clouds, land ice, sea 
ice, vegetation, and CO» (see Table 3). 

From Table 3, the equilibrium temperature change, AT.g=5.3°C and the temperature 
increment that would occur without feedbacks, ATo—0.6? C. The total feedbacks, У f; 
is related to the total gain, У) gi, by: 


1 | 
ARTE I oss ы (2.1) 


t 


and AT, is related to AT, by: 


AT 33° | 
AT, = 2 dod (2.2) 


eq = 
У fi 8.8 
The relative increment in solar energy required to yield the increment AT is: 
A Spo AT 4 x 0.6°K —3 
S ыы æ Vas) 


Note that the surface temperature of Earth is warmer by &33°K due primarily to the 
greenhouse gases. The change in solar energy required to obtain AT, is 


9.4 x 10? х 1367.54Wm ? = 12.8Wm ? 
that corresponds to a forcing at the tropopause level of 12.8 Wm" ?--5.7 or 2.2 Wm? 
This yields a climate sensitivity of 5.3/2.2 Wm"? = 2.4°C/Wm-?. 
Likewise, Hansen et al. (1984) estimate Y; f; at present to be 3 to 4, we take f ~3.6 


that yields a climate sensitivity of 1?C/Wm-?. If we assume in a world without ice, 
albedo effects are negligible, the temperature increase required to attain an ice-free state 
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FIGURE 3. The relation between equilibrium surface temperature (Teq), positive feedbacks 
(55 fi) and climate sensitivity. As a consequence of the increased positive feedbacks in colder 
climates, the same change in solar energy would have the greatest effect during glacial epochs 
and the least effect for an ice-free Earth. 


would be a minimum of 4.8°C (Barron & Washington 1984), the sum of feedbacks would 
be 2.1 and the climate sensitivity would be &0.6°C/Wm?. 

In order to go from the present climate to an ice-free (swamp-ocean) climate requires а 
forcing at the tropopause level of 8.6/Wm~? or an increase of 3.5% in the solar constant 
whereas a cooling from present to the climax of the last ice age at 18 ka requires only 
a forcing decrease of 2.2 Wm~? or a decrease in the solar constant of 9.4x 107? %o. 
Neither of those changes are obtainable by changes in solar irradiance alone (Table 2). 
However, significant solar modulation of climate change within the two steady states, 
glacial-inter-glacial, can be obtained. Neglecting the increased positive feedback during 
a Little Ice Age, which we estimate to be only 71076, a peak-to-peak change of 0.6°C 
(Table 2) explains most, if not all, of the Holocene temperature variation. The same 
change in solar irradiance, 2.4%0, would produce a peak-to-peak variation of 1.4°C at 
18 ka. Figure 3 illustrates the increased variability of climate during the ice ages and 
decreased variability at present and during an ice-free climate. 

The evidence that the recent warming event approaches or perhaps exceeds past events 
sets limits on natural variability during the Holocene but we must keep in mind that 
warming since the Little Ice Age is the combined effects of increased solar activity and 
greenhouse forcing. According to Hansen & Lacis (1990), greenhouse forcing to 1989 was 
2.04 Wm”?, to which must be added our estimated solar forcing of 0.58 Wm 7? (Table 2) 
or about 2.6 Wm"? total. This forcing at equilibrium, assuming the validity of the climate 
sensitivity of 1?C/Wm-?, would yield an increase of 2.6°C since the Maunder Minimum. 
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However, heat is transported to the mixed layer of the ocean and from there to the deep 
sea resulting in a long lag time. From their three-dimensional (3D) model, Hansen et 
al. (1984) estimated “the time for the global area-weighted mixed layer temperature 
to reach 63 percent of its equilibrium response in 124 years” (p. 156). The 124-year 
e-folding time is applicable for the high climate sensitivity of 1?C/Wm-?. If we divide 
the 0.6°C (Table 2) equilibrium temperature increase equally between warming after 
the Maunder Minimum and additional warming during the Modern Solar Maximum, we 
calculate that 25% of the current warming of 0.5? C is from solar forcing and 75% is the 
result of greenhouse gas forcing. The temperature rise following the Maunder Minimum 
would attain its equilibrium value of 0.3°C to which 0.125°C must be added from the 
warming due to the current solar maximum. The resultant 0.425? C would then constitute 
a little more than half of the 0.8? C warming since the Maunder Minimum. 

Our argument, thus far, has been based upon the high climate sensitivity of 1°C/Wm?. 
As a result of uncertainties in the feedbacks, the range of possible climate sensitivities lies 
between 0.5^C/Wm-? and 1.25?C/Wm-? (Hansen & Lacis 1990). The e-folding time 
is also a function of the feedbacks (У? f;). For the lowest value of climate sensitivity, 
0.5°C/Wm~?, the e-folding time is only about 40 years. This much shorter residence 
time slightly more than compensates for the lower greenhouse forcing climate sensitivity 
but the solar contribution falls to 1396 of the current global warming. The total tem- 
perature rise since the Maunder Minimum becomes 0.65?C, of which only 3396 is due to 
solar forcing. 

It is of interest to apply this analysis to pre-anthropogenic climate change. For this, 
we take predicted temperature change between the Medieval Solar Maximum (ad 1100- 
1250) and the immediately following Wolf Minimum (ad 1250-1350). We assume, as 
above, that the equilibrium warming and subsequent cooling would each be 0.3? C. The 
predicted temperature fall must be calculated from the e-folding time of 124 years. It 
would be &0.4°C. This global change is not great but would be observable in climatically 
sensitive regions using appropriate paleoclimate indicators. It would probably be below 
detection levels in other regions. If the climate sensitivity is 0.5°C/Wm~?, the tempera- 
ture decrease would only be &0.3°C but probably still observable. At high latitudes and 
at high elevations, the temperature might change by 1?C as compared to a 2°C warming 
during the current global warming event. There remains the possibility of greater cooling 
during prolonged periods of intense solar minima such as the Little Ice Ages (ad 1250- 
1850), for example. To completely account for Holocene temperature changes, however, 
it seems that all permissible solar and terrestrial factors must be maximized. 


3. Conclusions and summary remarks 


1. Pleistocene climate alternates between two steady states, glacial and interglacial, 
separated by transitional episodes. Climatic fluctuations around steady state are greater 
during glacial epochs than interglacial epochs by about a factor of three to four as 
observed in oxygen isotope stratigraphy. 

2. Regional temperature changes are much more variable than and should not be 
confused with global temperature change. Since climatic events like the Little Ice Age, 
Medieval Warm Epoch and Altithermal were originally defined based upon paleoclimatic 
and historical events from Northern Europe and the North Atlantic regions, there has 
been a tendency to exaggerate global temperature changes during the Holocene. 

3. Present global warming can serve as a frame of reference in evaluating Holocene 
global temperature changes from paleoclimate records within climatically sensitive re- 
gions such as high-latitude and high-altitude continental areas. Such comparison strongly 
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suggests that global temperature change did not exceed +0.5°C during the last 10 ka or 
+2.0°C during the last glacial period. 

4. The Sun is a low frequency harmonic oscillator with solar activity fundamental 
frequencies of 1/2120 yr, 1/212 yr, 1/87 yr and 1/11 yr plus harmonics. The 1/2120 
yr cycle includes two phases, episodes lasting 600 years during which the century-scale 
variation includes pronounced minima such as the Maunder Minimum separated by 1500- 
year episodes of high and more chaotic activity. 

5. The Sun in one of a family of Sun-like stars most of which have cycles of brightness 
similar to the Sun's. These Sun-like stars vary in brightness up to four times the variation 
in brightness of the Sun. Some of the Sun-like stars are in quiescent states similar to the 
Sun during the Maunder Minimum (Baliunas & Jastrow 1990; Lockwood et al. 1992). 

6. Lean et al. (1992) suggest that, during a solar minimum, the Sun decreases in 
brightness by 2.4%o relative to the present Sun similar to the average variation of the 
non-cycling stars (Table 2). 

T. Based on the 3D GCM analysis of Hansen et al. (1984), decreases in irradiance 
of 2.4% during solar minima relative to the present Sun would result in decrease of 
equilibrium global temperature by from 0.3°C to 0.6°C during the present interglacial 
and 0.9? C to 1.4°C during the last glacial maximum at 18 ka. 

8. Applying paleoclimatic restraints on global temperature and using the 3D GCM 
analysis of Hansen et al. (1984) with a climate sensitivity of 1°C/Wm~?, we estimate 
that global temperature increased by &0.3°C following the Maunder Minimum and more 
than 0.1°C of the 0.5? C of the current global warming since the late 19th Century would 
be the result of increased solar irradiance during the Modern Solar Maximum. Thus, 
about 5096 of warming since the Maunder Minimum would be the result of increase solar 
activity. If the lowest estimated climate sensitivity is assumed (0.52C/Wm-?), the solar 
contribution would be only about 33%. 

9. A climate sensitivity of 1°C/Wm~? requires that, if no more greenhouse gases were 
released to the atmosphere, global average temperature would continue to increase by 
2.6? C due &20% to solar activity and &80% to the effect of greenhouse gases. Such an 
increase would be unprecedented during the last 2 Ma. With inevitable further release 
of greenhouse gases, global climate would tend towards an ice-free globe held in check 
only be the large thermal inertia of the global oceans and polar ice. 

10. Assuming the validity of the astrophysics-solar physics analysis of Lean et al. 
(1992), changes in solar irradiance cannot result in a change from interglacial to glacial 
as suggested by Opik (1965). However, the century-scale and millennial changes in 
solar activity can account for much of the climate fluctuations around the glacial and 
interglacial steady states due to the existence of strong positive feedbacks. 

11. In order to completely explain Holocene global temperature changes by solar 
forcing, all terrestrial and solar factors must be optimized. Since this is unlikely, we 
conclude that solar forcing of pre-anthropogenic global climate change is a significant 
and perhaps dominant factor, but other processes must also be significant. 

12. Added note: The Vostok Antarctica core demonstrates relative climate stability 
during the present interglacial (Holocene) and also the penultimate (Eemian) interglacial 
(Dansgaard et al. 1993). The apparent instability during the Eemian suggested by the 
recent GRIP and GISP cores from Greenland (ibid.) appears to be caused by layer 
thinning and deformation rather than climate instability (e.g., Grootes et al. 1993). 
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L. L. Hood and J. L. Jirikowic 


Lunar and Planetary Laboratory, University of Arizona 


Tucson, Arizona 85721 USA 
ABSTRACT 


A 2200-2600 year quasi-periodicity is present in atmospheric A!*C records after re- 
moval of long-term trends due to the geomagnetic dipole amplitude variation. This period- 
icity consists of both a long-term variation of the mean and a superposed, approximately 
recurring pattern of century-scale variations. The strongest of these latter variations occur 
near maxima of the ~ 2400 year A!*C cycles. The residual record can be modeled to 
first order as an amplitude modulation of a century-scale periodic forcing function by a 
~ 2400 year periodic forcing function. During the last millennium, the largest century- 
scale variations (occurring near the most recent 2400 year A!*C maximum) are known 
to be mainly a consequence of the pronounced Maunder, Sporer, and Wolf solar activity 
minima, as verified by independent proxy solar activity records. Therefore, during this 
period, amplitude modulation has been occurring primarily in the sun and not in the ter- 
restrial radiocarbon system. It is therefore inferred that the ~ 2400 year forcing function 
is mainly solar although some secondary terrestrial feedback into the A1*C record is likely. 
This conclusion has implications for the predictability of future pronounced solar activity 
minima and for the interpretation of certain minor Holocene climatic variations. 


INTRODUCTION 


Deviations of atmospheric !*C versus time (A!*C) are produced in part by solar- 
induced changes in galactic cosmic ray flux which in turn modulate the radiocarbon pro- 
duction rate [Stuiver, 1961; Stuiver and Quay, 1980; Sonett, 1984; Damon, 1988; and 
references therein]. Directly dated tree ring A!*C records covering the last ~ 8000 years 
therefore represent one of the best available proxy measures of solar magnetic variability. 
However, an important problem in the interpretation of such records is the separation of 
solar variability contributions from other potential terrestrial sources of radiocarbon vari- 
ability. In this paper, empirical methods are used to investigate whether millennium-scale 
variations in the radiocarbon record are primarily solar or terrestrial in origin. 


ANALYSIS 


Figure 1 shows a high-precision A!*C record resulting from the 12th International Ra- 
diocarbon Conference [Stuiver and Kra, 1986]. A long-term trend is present that has been 
shown to be largely consistent with a modeled response to the ~ 11,000 year geomagnetic 
dipole moment amplitude variation as derived from archeomagnetic records [Sternberg and 
Damon, 1983; Damon, 1988]. Figure 2 shows the same record after detrending by removal 
of a least-squares-fitted cubic polynomial representing long-term changes in geomagnetic 
dipole moment intensity. Note that the residual radiocarbon variations are now plotted 
with an inverse scale so that decreases in ЛІС are upward and correspond to increases in 
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Figure 1. 12th International Radiocarbon Conference calibration record produced by com- 


bining high-precision records from a series of laboratories [Stuiver and Kra, 1986]. (Data 
courtesy of M. Stuiver and R. Kra) 
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Figure 2. Residual Conference А!*С record after detrending to remove the long-term ge- 
omagnetic variation and smoothing to minimize variations with periods less than approx- 
imately 100 years. Note that the vertical scale has been inverted so that A!*C maxima 
map solar activity minima. A least-squares-fitted sinusoid is superposed to indicate the 
approximate locations of maxima and minima of the ~ 2400 year quasi-cycle. 
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solar activity. 


The residual variations of Figure 2 exhibit a 2200-2600 year quasi-periodicity that has 
been recognized in power spectral analyses and visual examinations of various radiocarbon 
records [Houtermans, 1971; Denton and Karlén, 1973; Suess, 1980; Sonett and Finney, 
1989]. To indicate the approximate locations of maxima and minima of this quasi-cycle, a 
sinusoid with a ~ 2400 year period was least-squares-fitted to the residual record (fitted 
parameters: period = 2310 years; amplitude 4.08 per mil) and is superposed on the record 
in Figure 2. A geomagnetic origin for this quasi-periodicity in A!*C has been considered 
unlikely since there is no dipole moment amplitude variation near a period of 2400 years 
[Damon and Sonett, 1990]. However, there is a well-known climatic period near 2500 years 
[Mitchell, 1976; Pestiaux et al., 1987]. Several earlier investigators have reported that 
this ~ 2500 year climate variation has been nearly in phase with the inverted A!*C record 
during the last 8000 years [Bray, 1968; 1970; Denton and Karlén, 1973]. Figure 3 compares 
the residual A!*C record of Figure 2 to well-dated proxy northern hemispheric climate 
indicators including 6180 from the Camp Century Greenland ice core [Dansgaard et al., 
1984] and the Devon Island Canada ice core [Fisher, 1982] as well as to the temperature- 
sensitive Campito Mountain bristlecone pine tree ring width record of LaMarche [1974]. 
A tendency for climatic minima to be associated with A!4C maxima is evident. The last 
of these climatic minima was the Little Ice Age [Eddy, 1977]. 


There are two “end-member” candidate sources of the ~ 2400 year quasi-periodicity 
in atmospheric A1*C. First, the A!*C could be responding entirely to a terrestrial climate 
cycle of uncertain origin through a redistribution of 14C between the atmospheric and 
oceanic reservoirs. Second, the A!4C could be responding entirely to a solar quasi-cycle of 
the same period and the ~ 2500 year climate period could be independently driven by solar 
variability. In order to distinguish between these possibilities, it is helpful to consider the 
detailed characteristics of the radiocarbon record of Figure 1 together with independent 
proxy solar activity records covering the last 1000 years. 


Century-scale variations in the A!*C record are believed to be dominantly of solar 
origin as evidenced by their correlation with the independently verified Maunder and Spórer 
solar activity minima of the last millennium [Eddy, 1976; 1977; Stuiver and Quay, 1980; 
Stuiver and Braziunas, 1989]. As can be seen in Figure 2, the strongest century-scale 
variations tend to occur near successive maxima of the ~ 2400 year A!*C cycles. In 
addition, it has been shown that century-scale variations in the most recent ~ 2400 year 
cycle in A'*C are positively correlated with similar short-term variations in each of the two 
previous cycles [Hood and Jirikowic, 1990]. Thus, the quasi-periodicity consists of both 
a long-term variation of the mean and a superposed, approximately recurring pattern of 
century-scale variations. 

The behavior of the residual A!*C record can be interpreted in terms of amplitude 
modulation of a century-scale solar forcing function by a longer-term forcing function 
[Sonett, 1984]. Figure 4 shows a simplified example of amplitude modulation in which a 
200-year sinusoid is modulated by a 2400-year sinusoid. Pronounced short-term maxima 
(note the inverted scale) occur in the resulting time series (Figure 4c) at intervals of the 
longer period. If the 2400-year forcing function (represented in Figure 4b) were 
dominantly terrestrial in origin, then we would expect to be able to observe 
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Figure 3. Comparison between the residual Conference A!*C record and several well-dated 
proxy northern hemispheric climatic indicators for the last 8000 years (see the text). 
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Figure 4. Simple model of amplitude modulation of a 200 year sinusoid (a) by a 2400 year 
sinusoid (b). The product is shown in (c) and has characteristics that are qualitatively 
similar to those of the residual A!*C record. 
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Figure 5. À comparison between the last 1000 years of the residual Conference A1*C series 
and two independent proxy solar activity indicators (see the text). 
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the unmodulated solar signal (represented in Figure 4a) in independent proxy 
solar records covering the last 1000 years. However, the latter records show the 
occurrence of pronounced minima in the form of the Wolf, Spôrer, and Maunder solar 
activity minima. 

Figure 5 compares the residual A!*C record for the current millennium with the 
auroral frequency record compiled by Schove [1987] and with mean annual sunspot number 
since 1670 [Waldmeier, 1961; Eddy, 1976]. Within the recent interval covered by the 
sunspot record, both the Dalton and Maunder sunspot minima are approximately reflected 
in the residual inverted A!*C record. At earlier times, the auroral record remains positively 
correlated with the inverted A!*C time series (R = 0.55 at 20 years lag versus the smoothed 
auroral record for the entire series) and exhibits minima that can be identified with the 
Maunder, Spórer, and Wolf minima in the A!*C time series. A broad maximum centered 
on A.D. 1150, known as the medieval maximum, is also evident in both records. From 
the existence of pronounced minima in these independent proxy solar records 
near the last A!^C maximum, it is clear that amplitude modulation of century- 
scale solar variations has been occurring in the sun during the last 1000 years 
and not just in the terrestrial radiocarbon system. By inference, the earlier strong 
variations occurring at ~ 2400 year intervals are most probably produced mainly by solar 
changes although some secondary climatic feedback effects on A!*C can not be excluded. 
On this basis, the ~ 2400 year forcing function that is modulating the atmospheric A!*C 
record is suggested to be primarily solar. 


IMPLICATIONS 


Clearly, the construction of more complete quantitative models for A!*C production 
and exchange appropriate for millennium scale variations would be valuable and may 
elucidate further the relative roles of terrestrial and solar forcing mechanisms. However, 
the empirical evidence discussed here for a dominantly solar origin of both the century-scale 
and longer-term (~ 2400 year) residual variations in the Conference A1*C record does allow 
several provisional implications to be stated. The first of these relates to the predictability 
of future long-term activity changes. Although the Babcock-Leighton models for the Hale 
magnetic cycle (involving distortion of an initially poloidal field by differential rotation in 
the convection zone) have been developed, these models along with dynamo theories for 
the origin of the magnetic field itself are not yet sufficiently advanced to allow extensions 
to longer time scales [e.g., Gough, 1977]. From a purely empirical standpoint, however, 
it is reasonable to expect that the inferred quasi-cyclicity will persist in the future. In 
particular, it can be projected that large Maunder- and Spórer-type minima are not likely 
to occur again until ca. A.D. 4000 and that these will be preceded by a lengthy period of 
relatively high solar activity. In addition, the above analysis (indicating that the ~ 2400 
year A1*C periodicity is dominantly of solar origin) supports a possible solar origin for the 
Holocene climatic cycle of the same approximate period (Figure 3). 
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En termes clairs, nous sommes presque arrivés au pic de chaleur de 
notre période interglaciaire et le réchauffement climatique est dü au 
cycle de Milankovic et non à l'activité humaine. L'argument de la 
stratosphère froide est discutable, car les échanges de chaleur entre 
la stratosphére et le vide sont tels que l'augmentation de chaleur y 
est trés minime. De plus, il faut voir que cette activité de Sapiens est 
engendrée par le réchauffement! Vous n'imaginez pas un instant 
qu'une planéte ayant une gigantesque calotte de glaciation et un 
climat moyen sub-polaire pourrait abriter 12 milliards de Sapiens? Oü 
trouveraient-ils l'énergie pour se chauffer et cultiver de quoi manger? 
Rendez-vous dans 5000 ans! GT 
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Svensmark for ever! 


Il se pourrait que les tornades en augmentation qui proviennent de 
la vapeur d'eau en augmentation des océans tropicaux, soient en 
partie due à l'effet Svensmark. Rien ne dit que cet effet doit se 
produire d'abord dans les régions subpolaires. De plus, il est clair 
que nous sommes au pic de chaleur d'une des courbes de 
Milankovic, qui pourrait se terminer d'ici quelques milliers 
d'années. 


Un bon conseil: quittez la Floride, qui sera récupérée sous peu par 
les serpents, l'eau saumátre, les scorpions et les dolines, grace aux 
tornades et à la montée des eaux. Vous feriez bien de vous 
dépécher, mes mecs! 


Enfin je vous rappelle que la croissance de Sapiens en nombre est 
due au réchauffement climatique, et pas le contraire. 
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En termes clairs, nous sommes presque arrivés au pic de chaleur de 
notre période interglaciaire et le réchauffement climatique est dû au 
cycle de Milankovic et non à l’activité humaine. L'argument de la 
stratosphère froide est discutable, car les échanges de chaleur entre 
la stratosphere et le vide sont tels que l'augmentation de chaleur y 
est trés minime. De plus, il faut voir que cette activité de Sapiens est 
engendrée par le réchauffement! Vous n'imaginez pas un instant 
qu'une planéte ayant une gigantesque calotte de glaciation et un 
climat moyen sub-polaire pourrait abriter 12 milliards de Sapiens? Oü 
trouveraient-ils l'énergie pour se chauffer et cultiver de quoi manger? 
Rendez-vous dans 5000 ans! GT 


